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OO ' !• Introduction 
p ■ 

The counterintuitive Aharonov-Bohm (AB) effect fV\ represents one of the most widely 
^ ■ discussed issues of quantum physics. It predicts that a charged particle can be influenced 
■ by a magnetic field "even if the particle is nowhere in the region of nonzero field strength" 
|I2]|. This intriguing effect has alternatively been attributed to a nonlocal feature of quantum 
mechanics, to a hitherto unexpected direct physical meaning of an otherwise unphysical vector 
potential, to a topological cause and so on It is experimentally well verified JSHU, and 
its possible applications have attracted much interest recently, see e.g. E 13 H [H [H HOl. A 
generally accepted intuitive and physical understanding, however, seems to be still lacking. 

As Hamiltonian operator for an electron (without spin) in the presence of a magnetic 
field one takes 

1 2 



H 



P + eA(x) 



/2me (1) 



where P is the canonical momentum operator of the electron, me the electron mass and — e its 
charge. A(x) is a vector potential, which is not unique and is related to the magnetic field by 
B(x) = curl A (x). In classical physics such a vector potential is an auxiliary quantity without 
direct physical meaning. It may happen that a magnetic field vanishes in a region while the 
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Figure 1. Schematic sketch of the Aharonov-Bohm effect. Electron scattering by an infinite, 
perfectly shielded cylinder of radius R, with a homogeneous magnetic field Be^ inside, 
depends on the magnetic flux. 



allowed vector potentials do not. An example is given in fig. [T] with a constant magnetic field 
which vanishes outside an infinitely long cylinder of radius R. Since the integral ^ A(x) ■ dx 
over a circle around the cylinder yields the flux, A(x) cannot identically vanish outside the 
cylinder. 

A quick, heuristic, way to obtain the AB effect is to note that an electron may 
pass on either side of an impenetrable cylinder (see fig. [B, thereby picking up a phase 
{e/h) J A(x) • dx, with the integral taken over the electron's path. The resulting phase 
difference is (e/h)^, where $ is the magnetic flux in the cylinder, and this leads to field 
dependent effects for scattering. Related effects in magnetic fields without shielding had been 
pointed out earlier fTTl fT2l . 

A more detailed derivation [[II [19l describes electron scattering by eigenfunctions of the 
Hamiltonian ([Tj) with A(x) = {^/27rr)e^, where r = a/x^ + > R. Impenetrability 
is taken into account by choosing eigenfunctions with zero boundary conditions. Since the 
vector potential outside does not vanish, the eigenfunctions differ from those for the free 
Hamiltonian H^^^ = P^/2me with the same boundary conditions. This then gives the baffling 
result that electron scattering depends on the magnetic field inside the cylinder although the 
electron wavefunction cannot penetrate there. Ref. [3J explicitly calculates for a special case 
the momentum transfer to the electron from the wall of the excluded cylinder, compares it 
with the momentum transfer implied by the scattering cross section, and concludes that in 
general the force exerted by the cylinder surface is needed to satisfy Ehrenfest's theorem. It 
remains, however, physically somewhat puzzling why this force depends on the magnetic field 
inside the cylinder although the wavefunction vanishes on the surface and inside. Ref. |[T3l 
invokes a new, previously not considered, classical force to explain the AB effect. This was 
criticized in [T4] and rebutted in [15J. Ref. [[T6l draws a distinction between the Aharonov- 
Bohm phase shift and the Aharonov-Bohm effect and suggests that the Aharonov-Bohm phase 
shift is actually due to classical electromagnetic forces when relativistic effects are taken into 
account. In a recent experiment f,17J the time delay for an electron passing by a "macroscopic" 
solenoid has been investigated. 
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Some authors have noted that if one starts right away with a perfectly shielded cylinder - 
corresponding to quantum mechanics in a plane with a hole or in three-dimensional space with 
an excluded cylinder - the quantum dynamics for a particle outside is not uniquely determined 
(c/, e.g., [18J). In this point of view the choice of one particular quantum dynamics then 
appears somewhat ad hoc since it typically relies on information from the excluded region. 
This problem does not arise in approaches which consider finite, but increasingly high, 
barriers and then take a limit [fT9l [TSlI . In this case there are no interpretational problems. 
But it has also been shown [[T9l [T8l that as the barrier height is increased in the limit one 
arrives at exactly the same zero boundary conditions for the wavefunction as before. Thus the 
conceptual problem remains by which physical mechanism information about the inside field 
is transmitted to the outside. 

In this paper it is shown that in case of large, but finite, shielding the small wavefunction 
tails, which can tunnel into the cylinder and into the magnetic field, determine the force acting 
on the electron. In the limit of infinite shielding this force remains finite and flux dependent, 
it determines the otherwise under-defined dynamics outside and yields the AB effect. 

The paper is organized as follows. In section|2]we review the argument why the quantum 
dynamics can be viewed as ambiguous if one has a perfectly shielded cylinder right from 
the beginning. This is due to the fact that classically equivalent Hamiltonians can become 
physically inequivalent after transition to quantum mechanics. In section[3]expressions for the 
forces for finite and infinite shielding are given in terms of radial derivatives of the modulus of 
the wavefunction. These radial derivatives are calculated in section |4] and the resulting forces 
are determined. In section |5]it is indicated how the forces determine the dynamics and the AB 
effect. In the appendix expressions for the forces are explicitly derived which are valid for 
arbitrary geometries, not only cylinders. 



2. Indeterminacy of the quantization 

We first consider a free classical electron outside an impenetrable infinitely long cylinder of 
radius R. No assumption about a possible magnetic field inside the cylinder is made. By 
symmetry one can confine oneself to the x — y plane. A possible classical Hamiltonian 
function is H^^^ = P^/2mc, with reflections from the boundary of the excluded region 
implemented by the configuration space. The Hamiltonian, however, is not unique. Indeed, if 
^dum(x) is any vector function (a "dummy field") outside the cylinder, satisfying 

curlOdum(x) = 

there, otherwise arbitrary and not connected to any field inside the cylinder, then 

Hdum = [P + rJdum(x)]' /2me (2) 

equals the kinetic energy and also yieldfl i = 0. This is, up to a constant, the most general 
form of the Hamiltonian yielding x = and equaling the kinetic energy. Classically, all these 
Hamiltonians are physically equivalent. This general equivalence, however, is no longer true 

I This corresponds to the well known freedom to add to the Lagrangian a total time derivative. Cf., e.g., lITSl . 
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Figure 2. Wavefunction behavior near the boundary: finite shielding. Cylinder of 
radius R, wavefunction corresponding to incoming electrons of momentum hkex with kR = 
4.3 X 10~^. Two different magnetic fluxes: (a) $ = 0.2 h/e, and (b) $ = 0.4 h/e. Solid lines: 
|7/'^*'^''(r, ip) \ near the cylinder boundary for fixed polar angle Lp = I.Stt and for increasing 
barrier heights V as well as for the limit V co. Thick vertical line: cylinder boundary. 
Broken lines: slope at this point. For the two fluxes the slopes differ significantly, but differ 
little for a fixed flux and different barrier heights V . 




Figure 3. Infinite shielding with dummy field. kR as in fig. |2] (a) SoHd lines: |?/'['^ (r, 
near the boundary as a function of r for fixed polar angle if = l.Svr and for two different values 
of K which labels inequivalent dummy fields A^j'^^^-'^^^ from (|5]l. The slope at r = i? (dashed, 
dot-dashed lines) depends on k. (b) The slope at ?■ = i? plotted as a function of k for two polar 
angles tp^ — n and ~ I.Stt. The value of k is uniquely determined if the slope is known 
for two different angles Lp. 



after transition to the quantum theory (with zero boundary conditions on the cylinder) {cf, 
e.g., [[I8l|20l). For a given f2dum(x) with curl r2dum(x) = outside the cylinder, we define k 
as 

K = non — integer part of <j) r2dum(x) ■ dx//i (3) 

so that < K < 1 . Then one can define the continuous and differentiable function of modulus 
1 

■dx'l (4) 



A(x).exp{l/ 



^dum(x ) Gipi 
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where xq can be chosen as xq = (-R, 0); r, r' > R. It is easy to check that A(x)iydumA(x)^ 



^duL where 



rV(«:) 
dum 



/2me, flil(x) = -e, (5) 



with < K < 1, and the zero boundary conditions are preserved. The Hamiltonian operators 
-^dum physically inequivalent for < k < 1, and without further input information 
it is unclear which one to choose. Thus there is an ambiguity when one starts right away 
with infinite shielding, and this ambiguity can be fixed by imposing an additional boundary 
condition which arises in a natural way when one takes into account tunneling and associated 
forces, as seen in the following. 



3. The forces for finite and infinite shielding 



Finite shielding of the cylinder can be modeled by a potential V(x). Just for simplicity we 
consider V(x) = Vo9(-R — r), with Vq — oo later on, and a homogeneous magnetic field 
inside the cylinder with vector potential 



A$(x) 



2^ 



^e(i?-r) + le(r-i?) 



X 



(6) 



where $ is the magnetic flux, but also other magnetic fields can be considered. The 
Hamiltonian then is 

P2 + eA. 



L$ (X 



/2me + V (x) 



(7) 



We consider a normalized (planar) wavefunction, ijji ' ' , corresponding for t — oo to an 
incoming free particle which is then scattered. Due to tunneling small tails enter the cylinder. 
The total force on the electron is 



(8) 



Note that only the tails contribute. As the tails, the Lorentz force goes to when V oo. 



For a step potential the first part becomes Vq ^ Rdip l-ipl ' {R, V^) pe^. Expanding tjjj: ' in 
terms of eigenfunctions of His>,v and using eqs. (35) - (40) of [191 one can show explicitly for 
large Vo that V^^^\4'''^\r, if)\ = hdr\4'''^\R, </')l/(2m)V2 + 0(K,"')- Thus 



2mp 



2tt 



Rdip 



d_ 

dr 



1 2 



e. + 0(V)- 



(9) 



This result is also true for more general V and B, e.g. B(r) = 0, r > R/2, as well as for the 
torus and other domains (with a surface integral and normal derivative). This general case is 
treated in the appendix. 



For infinite shielding and fixed arbitrarily chosen dummy field Q 



d/dt {ip^'^^l'me'vliljl'^') for the force and eventually obtain 



dum we use F, 



2m„ 



2n 



Rdip 



d_ 

dr 



(10) 



Again this carries over to other domains as shown in the appendix. 
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The as yet underdetermined dummy field can be made unique by requiring that the 
K-dependent force (fTOl) agrees with dH) in the limit V ^ oo (obtained with the field actually 
contained in the cylinder). It will now be shown that this requirement uniquely determines k, 
and thus also the dynamics of the idealized case. 



4. The radial derivative and resulting forces 

It suffices to consider electron beams of definite incoming kinetic momentum, hkex say. 
The corresponding scattering solution is a stationary state, denoted by V'fc*'^'* and V'^'^'*, 
respectively. To calculate these we note that around the backward direction, ip = -k, the 
incident wave should behave as an eigenfunction of Pkin- This implies that the incident wave 
should behave as 

gik.x-*e${^-^)/h kr^l, \^-'K\<e (11) 

(k) 

and similarly for . The asymptotics in Ref. [JJ, based on the probability current, is the 



sami _ 

The wavefunction %lj'f''^\r,ip) can be calculated analytically as a series and evaluated 
numerically. Note that the wavefunction is gauge dependent while its absolute value is not. 
Although the tails go to zero inside and on the cylinder when V tends to infinity, the behavior 
of the wavefunction in the vicinity of the cylinder depends sensitively on the value of 
as seen in figs. |2] (a) and |2] (b). There we consider two different magnetic fields inside 
the cylinder, with fluxes $i and $2, and increasing barrier heights V . For each flux and 
barrier height ilj'f'''^\r, (p) is calculated and its absolute value is plotted as a function of r for 
y9 = 1.3 TT. With increasing barrier height V the wavefunctions are indeed seen to converge to 
zero inside the cylinder (r < R). The form of the tails depends on the specific B and V, while 
the limiting slope depends only on $, or rather on 

a = non — integer part of e$//i. (12) 

With infinite shielding and a dummy field rij^^^^^, the scattering solutions ilj[^\r,(p) 
vanish on the cylinder (r = R), but the rate with which is approached when r ^ R 
depends on k. To see this we calculate ipl^^ (r, ip) by expressing it as a linear combination 
of eigenfunctions of H^J^ of fixed canonical angular momentum, with unknown coefficients. 
Using the asymptotics of Bessel and Hankel functions and (fTT)) one can determine the 
coefficients and obtains 

00 



J\n+K\ikR) 



J\n+K\{kr) - ^— -iJ|„+^|(/cr) 



(13) 



§ One can also start with an incident plane wave |k) ^e*'' and use the Lippmann-Schwinger equation to obtain 
the corresponding scattering solution. We have determined the necessary Green's functions and have shown that 
one obtains the same expression as with (fTTl l. This generalizes a result for a magnetic string [21 J and will be 
published elsewhere. 
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The rate is given by the tangent slope of l?/'^"'' (r, | at the cylinder, i.e. by the radial derivative 



-R, ip)\ which can be calculated from (fT3l) . Since ip^^^ = on the boundary of 



the cylinder one easily sees that 



dr 



if 



R. 



(14) 



With this identity and using the Wronskian for Bessel and Hankel functions, one obtains from 
(fni) 



dr 



2 



E 



Jn(Lp+-K) 



We note that this derivative is invariant under the unitary transformation A(x) in dH). 



(15) 



In fig. [3] (a) \ip'^^\r, v^)] is plotted as a function of r for two values of n and for fixed 
Lp = I.Stt and it is seen that the slope at r = i? depends on k. In fig. [3](b) the slope at 
r = i? is plotted as a function of k for two values of ip (solid curve: ipi = n, dashed curve: 
(p2 = I.Stt) and it is seen that the values of the slope at the two different angles determine k 
uniquely. Moreover, this slope agrees with the limiting slope for a magnetic flux $ if k equals 
the non-integer part of e^/h (i.e. if k = a), and then the corresponding forces on the electron 
beam are equal, by eqs. dH) and (flOl) : conversely the condition k = ais also necessary for this. 
Thus requiring the equality of the forces determines the previously underdetermined dummy 



dum' 



field n 

The force on the electron beam can easily be calculated for V — > oo by means of 
(fT5l) (with K = a) as a rapidly converging sum involving Hankel functions. There is both 
a component in the backward direction as well as a perpendicular component. The former 
is repulsive and invariant under the replacement a 1 — a. The latter reverses sign under 
a ^ I — a and under charge reversal, vanishes for a = 0, 1/2, 1, and for small $ points in 
the same direction a Lorentz force would do for an electron inside the cylinder. For kR -C 1 
one obtains for the force per unit cylinder length on an electron beam of incoming momentum 
hke^ and density p 



p(<i>,y— ►oo) 



P 



h^k ( -2sin2(7re<l>//i) 



OiR" 



(16) 



sin(27re<l>//i) 

where a' = min(2a, 4 — 4a). The higher-order terms in R contain the reflecting force by the 
cylinder and they vanish for R ^ (magnetic string). For kR ~ 1, F2 is several orders of 
magnitude smaller than Fi since reflection by the cylinder dominates. Subtracting this gives 
about the same order of magnitude, but the remainder is overall much smaller than in (fT6l) . 



5. Forces and the AB effect 



These results yield a physical explanation of the AB effect as follows. For the idealized case 
of an impenetrable cylinder, the quantum (but not the classical) dynamics in the outside region 
is underdetermined since it contains a largely arbitrary "dummy" field. Although tempting, 
there is no a priori reason to relate this field to an (in principle unknown) magnetic field 
inside the cylinder. This leads to physically inequivalent Hamiltonians h]^^^, < k < 1, 
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and a k dependent force. To derive an additional boundary condition which removes the 
indeterminacy we consider high, but finite, barriers. Then, by tunneUng, minute tails of the 
scattering wavefunction can enter the cylinder. The shielding and a magnetic field inside 
the cylinder exert a force on these tails and on the electron. By the combined influence of 
magnetic field and increasing shielding on the form of the wavefunction this force remains 
finite and flux dependent even when the shielding goes to infinity and it can be expressed by 
the radial wavefunction derivative at the cylinder. Requiring that the force for the directly 
tackled idealized case (with infinite shielding from the outset) be equal to this limiting force 
fixes the former's as yet underdetermined dynamics {i.e. k). Alternatively, the limit slope of 
the scattering solution at the boundary can be considered as an additional boundary condition 
for the ideal case which also removes the indeterminacy. The dummy field determined in this 
way is, up to a factor e, just the vector potential customarily used right away in the discussion 
of the AB effect. 

6. Summary 

It has been shown that the AB effect for cylinders arises quite naturally when one considers 
tunneling and the force exerted on the small wavefunction tails of the electron inside the 
cylinder. Although the Lorentz force vanishes when the tails go to zero in the limit of infinite 
shielding, the flux dependence of the remaining force persists and precisely yields the AB 
effect. It has also been shown that the limit slope of the scattering solution at the boundary 
can be considered as an additional boundary condition for the ideal case which also removes 
the indeterminacy of the quantization procedure. The same results are expected to carry over 
to the torus and other domains. 

Appendix A. General case: Forces in a magnetic and scalar field 

Appendix A. 1. Magnetic and large, but finite, scalar potential 

We first consider a general time-independent magnetic field B(x) in a region G, with a vector 
potential A(x) vanishing at infinity, and a scalar potential V{x.) which is nonzero in the same 
region G and vanishes outside. Later we will remark on the more general case that the scalar 
potential may also vanish on parts of the interior of G. 

The force, Ff^\ on an electron of charge — e at time t is then the sum of the scalar and 
Lorentz force. 



where V't denotes the wavefunction under the time development with the Hamiltonian Hy, 
of an electron coming in from infinity. The Hamiltonian is 




(#\(-VV--evxB)#)) 



(A.1) 



Hv^(P + eAf/2m + V . 



(A.2) 
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We will investigate the behavior of the force in (lA.ll) for large scalar potential V and will 
show that it can be expressed as an integral over the surface, dG, of the region G, 

2 





[ do 




2m 


JdG 


dn 



(A.3) 



+ terms — + for V 



oo 



Note that the dependence of this expression on the magnetic field comes through the time 
development via the Hamiltonian. In two dimensions and if the region G is a circle one 
arrives at Q if one lets ip^^^ go to a stationary solution of Hy. 

To prove (IA.3I) we first consider the scalar part in (|A.1I) and show that 



2m jQ 
+ terms 

To show this we use partial integration to write 



+ CO. 



for y ^ oo 



(A.4) 



(V) 



vi^^;^^) = / <fxv{^) V 



(A.5) 



G 



V can now be expressed by the Hamiltonian in (IA.2|) . For the latter we have 



^((P + eA)#,"-'. 

1 



(V) 



((P + eA)#\ (P + eA)^rOiRAc 



2m 



(A.6) 



+ -((P + eAW,' 



[P + eA)#))G+(V^. 



G 



where the indices denote integration over the respective regions. The term (V't^^ Hyipl^^) is 
independent of t and therefore equals the incoming kinetic energy. The first term after the last 
equality sign converges for V ^ oo to the corresponding expression with Dirichlet boundary 
conditions on G, by standard results [[22l[T8l . e.g. if A is bounded and if the initial incoming 
wavefunction is in the domain of the Hamiltonian. This limit also equals the incoming kinetic 
energy. Since all terms are non-negative, the remaining terms in (IA.6I) have to go to zero 
when V ^ oo. Hence for ^ oo 

((P + eA)#\ (P + eA)#V -0 
With this one obtains for bounded A and with the inequality 



(A.7) 
(A.8) 
(A.9) 



[P + eA)^, 



iv) 



P^, 



iv) 



eA ijj- 



iv) 



\G 



that also 



P^. 



\G 



for y — i> oo 



(A. 10) 
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This and Schwartz's inequality then imply that 

(Wf^ Hv^Ij^P)g ^ for 1/ ^ oo (A.ll) 
if the initial wavefunction is in the domain of Hy- Inserting 

V = Hv-{^'^ + eKfjlm 

into (|A.5I) one obtains from (lA.llI) 



+ c.c. 

+ terms for ^ oo . (A. 12) 

Using (lA.lOl) and Schwartz's inequality this yields the claim in (|A.4I) . 

We will now investigate the first term on the right hand side of (IA.4I) and will show that 
it can be expressed as a surface integral, i.e. for the i— th component 
1 



— / d^x d.^^Pv^'^P + c.c. 
2m .]q 



2m 



di^P^V^P ■ do + c.c. 



- ■ / do V^P^ ■ Vtpp^ } . (A. 13) 



IdG 

To prove this we use partial integration, i.e. Gauss's theorem in the form 



/ fd,g = - [ {d,f)g + [ e,-dofg. (A. 14) 

Jg Jg JdG 



Then 



d'^xd^ijPdjdjijP = - / d'^xdjdiijPdj4 



G JG 



+ I ej- do diijPdjijp^ (A. 15) 



dG 



and 



d'^x diipPdjdjtpP = - d'^x djditpPdjtpP 



G JG 



+ I eydod.ijPd.ijP 



dG 



(A.16) 



Applying (|A. 141) to the first term on the r.h.s. of (IA.16I) yields 



d'^x ditpPdjdjtpP = / d'^x djtpPdidj^pP 



G JG 



— Gj ■ do djtpPdjtfjP 

JdG 



+ I ej- do ditpPdj^''^ 



dG 



(A.17) 
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Adding eqs. (IA.15I) and (IA.17I ). the first terms cancel and one obtains 



c.c. 



G 



JdG JdG 

— I ■ do djipl^^ djipl^^ 

JdG 

(A.18) 

which gives (|A.13I) . 

For large V (IA.13I) can be further simplified. Indeed, one has i/jI^^ — on dG for 
V ^ CO and hence 

V^^^ — - — n (n = normal) 

on 

rj , {V) 

a.#)= e..V#)-e..n^ (A.19) 

on 

for — > oo. Insertion into (IA.13I) gives 
1 



— / d^'x dii^Y "P'^t + c.c. 
2m 



G 



= 7;— ■ ndo— — (A.20) 

2m Jqq on on 

plus terms going to for V ^ 00. 
For the Lorentz force one has 

and one finds from (|A.10I) and Schwartz's inequality together with the boundedness of B that 
it decreases to when V ^ 00 . This, together with eqs. (|A.20I) and (IA.4I) yields (IA.3I) . 

These results can be generalized to the situation when the potential V vanishes on a part 
of the interior of the region G, e.g. one could have a magnetic field in a cylinder of radius R 
while the shielding would be only in a cylindrical ring Ri < r < R. Then for particles coming 
in from infinity one obtains the same results as above since the incoming wavefunction will 
less and less penetrate through the ring. This can be made mathematically precise by a similar 
reasoning as above. 



Appendix A.2. Forces in the case of infinite shielding 

If one starts right away with infinite shielding the particle motion takes place completely 
outside the region G and there is no regular expression for the force acting at the boundary. 
We therefore use d/dt (mv) as definition for the force. Moreover, the Hamiltonian is not 
uniquely determined, as pointed out before, and we arbitrarily pick for it a particular -ffdum 
with a dummy field ildum, analogously to Q. Then 

mv = P + ridum- (A.21) 
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Now let ipt be a normalizable wavefunction coming in from infinity (and being in the domain 
of i/dum)- We will show that then 



{ipt, mw^t) = 7;— / do 



dG 



(A.22) 



dn 

Note again that the dependence of this expression on the dummy field Odum comes through 
the time development via the Hamiltonian. In two dimensions and if the region G is a circle 
one arrives at (flOl) if one lets ipt go to a stationary solution of Hdum- 
To show (IA.22I ) we write the left hand side as 

d ITfi 

— (V't, mvipt) = — {(i/dum^t, vV^t) - {ipt, vifdumV't)} (A.23) 

At this point it is crucial that in the first term i^dum cannot be moved over to the other side 
by hermiticity since then one would be led to [//dum, v] = x = and thus there would be no 
force. The underlying mathematical reason why this is not allowed is that vipt need not lie in 
the domain of Hdum- However, one can move Odum ■ P over by partial integration since the 
boundary terms vanish. Thus 

d iTTi f r 

-{il)umwil)t) = - — - / d'^x iP'^^ptVi/jt 

+ ^(2 0dum-P + riLm)vV't 

-V^vi/dum^J . (A.24) 



dt 



Inserting = ■iptfi'^V'^ {-vipt) — iptfi'^V^ {'vipt) gives 
at Zmrij]iid\Q 

Note that the last commutator vanishes. Inserting (IA.21I) the ridum terms cancel, by partial 
integration as in (IA.14I ). since the boundary terms vanish. Thus one obtains for the i— th 
component of the force in (IA.23I) 



dt 2m ./Kd\G 



2m 



d^'xVtPfVd.tPt- / do -V^PAi^t 

M''\G JdG 



+ / d'xV^Pt■V^,^Pt+ / do-^,V5,V't 

JjRd\G JdG 

where the last equality results from partial integration. The first and third term cancel, while 
the last one is zero since 'tp vanishes on the boundary. This yields 

^{iJt, mwi^t) = 77- /" do ■ V^VV^t . (A.25) 
dt 2m Jqg 

Since Vipt is perpendicular to the boundary dG this gives (IA.22I) . 
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